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The fully implicit ablation and thermal response code simulates pyrolysis and ablation of thermal protection 
materials and systems. The governing equations, which include energy conservation, a three-component 
decomposition model, and a surface energy balance, are solved with a moving grid. This work describes new modeling 
capabilities that are added to a special version of code. These capabilities include a time-dependent pyrolysis gas flow 
momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite- 
rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with 
these new additions. Two groups of parametric studies of the phenolic impregnated carbon ablator are performed. In 
the first group, an Orion flight environment for a proposed lunar-return trajectory is considered. In the second group, 
various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand 
the effect of pyrolysis gas momentum transfer on material in-depth thermal responses with finite-rate, equilibrium, or 
frozen homogeneous gas chemistry. Results indicate that the presence of chemical nonequilibrium pyrolysis gas flow 
does not significantly alter the in-depth thermal response performance predicted using the chemical equilibrium gas 
model. 
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Nomenclature 

permeability slip parameter, Pa 

preexponential constant, s -1 

m/p e u e C M , dimensionless mass-blowing rate 

Forchheimer coefficient 

Stanton numbers for heat and mass transfer 

mass fractions of species i 

specific heat, J/kg ■ K 

diffusion coefficient. m 2 /s 

diameter of pore, m 

total internal energy, J /kg 

internal energy, J/kg 

activation energy, J /kmol 

total enthalpy, J /kg 

enthalpy, J /kg 

recovery enthalpy, J /kg 

wall enthalpy, J/kg 

permeability, m 2 

thermal conductivity, W/m ■ K 

permeability in the limit of continuum flow, m 2 

backward reaction rate 

forward reaction rate 

molecular mass, g/mol 

total number of homogeneous reactions 

mass flux, kg/m 2 • s 

total number of chemical species 

pressure. Pa 

conductive heat flux at surface, W/m 2 
radiative heat flux at surface, W/m 2 
universal gas constant, J/kmol • K 
gas constant, J /kg ■ K 
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= mass generation rate of gas species i due to resin 
decomposition, kg/m 3 • s 

= mass generation rate of gas species i due to 
homogeneous reactions, kg/m 3 • s 
= pore diameter-based Reynolds number 
= temperature, K 
= time, s 

= local gas velocity, m/s 
= local grid velocity, m/s 
= surface absorptance 
= volume fraction of resin 
= porosity 
= surface emissivity 
= blowing reduction parameter 
= viscosity, N • s/m 2 
= total density, kg/m 3 

= original density of pyrolysis gas component, kg/m 3 
= residual density of pyrolysis gas component, kg/m 3 
= Stefan-Boltzmann constant, W/m 2 ■ K 4 
= mass fraction of virgin material, defined in Eq. (11) 
= decomposition reaction order 


Subscripts 
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char 

boundary-layer edge 
pyrolysis gas 

density component (A, B, and C) or gas species 

surface species 

virgin 

wall 
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I. Introduction 

T HE authors have developed a family of programs for analysis of 
ablative thermal protection system (TPS) materials. The fully 
implicit ablation and thermal response code (FIAT) [1], the two- 
dimensional implicit thermal response and ablation code (TITAN) 
[2], and the three-dimensional finite-volume alternating direction 
implicit ablation and thermal response code (3dFIAT) [3] simulate 
the internal heat conduction, in-depth thermal decomposition, quasi- 
steady pyrolysis gas flow, and surface ablation of TPS materials in 
one, two, and three dimensions, respectively. FIAT is widely used by 
NASA and the industry as the one-dimensional analysis and sizing 
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Table 1 List of pyrolysis gas chemical reactions and rates 


Reaction number 
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tool for spacecraft TPS materials. TITAN can analyze problems with 
two-dimensional or axisymmetric geometry. In some cases, a two- 
dimensional analysis is inadequate, and a three-dimensional ablation 
code is required to perform a high-fidelity simulation. The 3dFIAT 
program can analyze the thermal response of the entire heat shield of 
a space vehicle. The prediction of ablative heat shield response for a 
spacecraft entering the atmosphere with an angle of attack is such 
a case. 

In this family of programs, decomposing materials are assumed to 
behave as three independently pyrolyzing components. Once the 
pyrolysis gas is formed, the gas flow is assumed to reach steady state 
instantaneously by ignoring the pressure and viscous resistance 
effects. These assumptions are generally good if the char depth is 
relatively thin. In such cases, there is no need to solve the momentum 
and species conservation equations for pyrolysis gas flow. 

The nonequilibrium chemistry of pyrolysis gas inside the char 
layer was first studied by April et al. in 1971 [4], Since then, there 
have been no significant efforts made in this area to model the finite- 
rate homogeneous chemical reactions of pyrolysis gas. Recently, a 
study conducted by Ayasoufi et al. attempted to examine the effect of 
nonequilibrium pyrolysis gas chemistry on the performance of 
charring ablators [5], in which the pyrolysis gas chemistry was based 
on the work of April et al. [4] . However, in the work by Ayasoufi et al. 
[5], the gas flow momentum conservation was not considered, and 
thus, the mass flow rate and kinetic energy of the gas could not be 
accurately predicted. The effect of permeability on gas flow inside a 
charring ablator was extensively studied by Ahn et al. [6]. They 
solved the gas momentum equation within a carbon-phenolic heat 
shield of the Pioneer Venus probes for pyrolysis gas in chemical 
equilibrium. Recently, the non-Darcy behavior of pyrolysis gas in a 
thermal protection system also was examined by Martin and 
Boyd [7], 

Given the large surface area within porous TPS materials, one 
would think that both homogeneous (gas) and heterogeneous 
(surface) reactions may be important for how quickly pyrolysis gas 
approaches a state of chemical equilibrium. Heterogeneous reactions 
are important for additional complex phenomena, such as coking, 
which was observed in the Apollo heat shield [8]. However, the effect 
of heterogeneous reactions within pyrolyzing heat shields has been 
mostly neglected. 

Darcy’s law has limited applicability, because it only accounts for 
the viscous resistance of the flow at a low speed. If the flow speed is 
not sufficiently low, then inertial contributions to the flow resistance 
may become noticeable. The Darcy-Forchheimer equation [9] is 
widely used to account for the inertial effects in porous media. In this 
paper, we enhanced the FIAT code by including the time-dependent 
pyrolysis gas momentum equation with the Darcy-Forchheimer 
terms within a porous char layer to obtain the flow velocity. Addi- 
tionally, the multispecies mass conservation equations for pyrolysis 
gas are solved to simulate the nonequilibrium chemistry that may 
occur inside a char layer. The three-component decomposition 
model, which is based on the thermogravimetric analysis (TGA) data, 
is used to estimate the instantaneous density of resin. For 
nonequilibrium pyrolysis gas modeling, the homogeneous chemical 
reactions and associated rates similar to those developed by April 
et al. [4] are considered. 

The charring thermal protection material considered in this work is 
the phenolic impregnated carbon ablator (PICA) [10], which is a low- 


density ablator that was used as the heat-shield material of the 
Stardust sample return capsule [11]. The purpose of this paper is two- 
fold: first, to describe the details of the formulation considered in the 
enhanced version of FIAT; and second, to perform parametric studies 
to understand how the pyrolysis gas momentum transfer affects PICA 
material in-depth thermal response with various gas chemistry 
models to guide the direction of future work. The aerothermal 
environments used for the parametric studies are those relevant to the 
Orion Lunar return mission, including both flight and arcjet 
conditions. The predictions are presented and compared with those 
based on the PICA/FIAT material response model developed by the 
Orion Thermal Protection System Advanced Development Project 
(TPS ADP) [12], 


II. Governing Equations 

Under the assumption of local thermal equilibrium, the total 
energy conservation equation in conservative form for the combined 
solid-gas system inside a charring ablator is defined as 

d 

— (pc p T + eE g ) 4- V ■ ( eu g H g ) + V • (kWT) 

= \ ■ V(pc p T + eE g ) (1) 

where 

V | V | 

E g = J2p‘ e ‘ + ^Ps u l H s = J2 p ' hi + 2 PgUlg 
1 = 1 1 = 1 

/ N s \ N s 

P = ( Pi R i ) T = Ps R J Pg = Y.P- 

M= 1 ' i=l 



Fig. 1 Aerothermal environment for group I. 
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Time (s) 


Fig. 2 Surface and in-depth temperatures predicted by baseline FIAT 
and by FIAT with momentum equation. 


The individual terms in Eq. (1) are interpreted as follows: rate of 
storage of total internal energy, net rate of total enthalpy converted by 
gas flow, net rate of thermal conductive heat flux, and convection rate 
of total internal energy due to coordinate system movement. The 
conductivity of pyrolysis gas is much smaller than that of the solid, 
and thus, the conductive heat flux through the gas is ignored. If the 
internal energy, p,-e ; , and the kinetic energy, ^p g u~, of gas flow are 
ignored, then Eq. (1) reduces to the same form used in the previous 
work [1], 

Darcy’s law is applicable to low-speed flows, such as the modeling 
of underground flows. Typically, this means that the Reynolds 
number, based on the average pore diameter, has an order of 
magnitude near unity or less: 



f* g 



0 200 400 600 800 1000 1200 

Time(s) 


Fig. 3 Pyrolysis gas mass injection rate predicted by baseline FIAT and 
by FIAT with momentum equation. 

Darcy’s law only accounts for the viscous resistance of the flow. As 
the flow velocity is increased, inertial terms also should be 
considered. The Darcy-Forchheimer equation is widely used to 
account for the inertial effects. The transient Darcy-Forchheimer 
equation adopted in this work to predict the pyrolysis gas flow 
velocity is written as 

dep„ u„ 

+ v ’ ( e Pg a g^g + e P) 

= - e Y U s- e ^Pg\ u g\ u g + ' f - S 7 e Pg ll g ( 2 ) 

where the first term on the right-hand side is the Darcy term and the 
second term is the Forchheimer term, and K and C F are the 
permeability and the Forchheimer coefficient, respectively. The third 
term accounts for grid motion in the numerical solution. The left- 
hand side of Eq. (2) is the unsteady form of the Euler equation with 
porosity. 



a) 


b) 

Fig. 4 Temperature and pressure profiles at times equal to a) 80 s and b) 515 s. 
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Fig. 5 Gas and solid density distributions at times equal to a) 80 s and b) 515 s. 


Most TPS materials are nonisotropic, and therefore, the quantities 
k and K are actually tensors. However, in one dimension, we may use 
the scalar quantities as written inEqs. (1) and (2). More general forms 
that include tensor material properties and Brinkman's term may be 
required for solutions in two or three dimensions. 

The permeability and the Forchheimer coefficient for both virgin 
material and char must be defined to perform a simulation. Some 
models and data for permeability and the Forchheimer coefficient are 
available. If the pressure is not too high, the effective permeability 
may be expressed in the Klinkenberg form as 



where K 0 is the gas permeability in the limit of continuum flow, and b 
is a permeability slip parameter that accounts for rarefied flow effects. 


Values of K 0 and b for the through-the-thickness direction of PICA 
are found in the work of Marschall and Cox [13]. The formula for the 
Forchheimer coefficient used by Ward [14] and Ahn et al. [6] is 
adopted. 

The pyrolysis gas is assumed to have N s components. The mass 
conservation for each component is 

+ v ■ epiU g + ep g D, V = e R. + R dl + v • Ve p t (4) 
and the global mass conservation is 

•j N„ 

+ v • (ep g u g ) = ^ fj di + V • Vepg (5) 



Fig. 6 Kinetic energy density and mass flow rate distributions at times equal to a) 80 s and 515 s. 
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Here, R di is the mass generation rate of gas species i due to resin 
decomposition, and Rj is the gas species i production rate due to 
homogeneous gas reactions. 

For pyrolyzing TPS materials, a standard three-component 
decomposition model is used. The instantaneous local density of the 
composite is given by 

P = r (pa + Pb ) + (1 _ Hpc (6) 


where the parameter T is the volume fraction of resin and is an input 
quantity. The three components decompose independently following 
the Arrhenius type relation 


^=-B s 

dt 


^W) p ° '^7rT +v ' vpi (7) 


where p ti is the residual or terminal density of component i, and p oi is 
the original density of component i. B a , E a , and V are determined 
using curve fitting of thermogravimetric analysis data. The mass 
generation rate of each gas phase chemical species, R Ai , formed 
during the resin decomposition can be estimated by multiplying the 
mass production rate computed in Eq. (7) by the mass fraction of each 
individual chemical species. If the resin decomposition process is 
under the chemical equilibrium condition, the mass fractions of 
chemical species can be computed using a chemical equilibrium 
solver, such as the aerotherm chemical equilibrium (ACE) [15] or 
multicomponent ablation thermochemistry (MAT) [16] codes. If the 
decomposition process is not in chemical equilibrium, details of the 
heterogeneous chemical reactions and their associated rates involved 
in the decomposition process are required to estimate the mass 
fraction of each chemical species. For the parametric studies 
presented in this work, the resin decomposition process is assumed to 





Fig. 7 Surface and in-depth temperature predicted by a) 119-species and nine-species chemical equilibrium models; b) finite-rate, equilibrium, and 
froSurface and in-depth temperatures predicted by zen gas chemistry models; and c) baseline FIAT and by enhanced FIAT with 9-species chemistry. 
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Fig. 8 Mass fraction profiles predicted by finite-rate and chemical equilibrium chemistry models at times equal to a) 80 s and b) 515 s. 
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be under the chemical equilibrium condition. The mass fractions of 
chemical species formed during the resin decomposition are 
predetermined and then arranged in tabular form as one of the FIAT 
input files. 

For a general homogeneous chemical reaction, the mass 
production rate R t of gas species i has the following form: 

m r N s N s 

Ri = Mi J2(Pij ~ r u ) k fj n ( c >') p ' - k n II ( c ;) r ' ( g ) 

7=1 L i=l 1=1 

where p, ; - and r,-,- denote the stoichiometric coefficients of species i as 
a product and a reactant of reaction j; m is the total number of 
reactions and N s is the total number of species; is the molecular 
mass; Ktj and K r j are the forward and backward reaction rates; c ; is 
the species concentration; and p[ and r\ are the orders of reaction 
with respect to species i. 

A complete set of homogeneous pyrolysis gas chemistry for PICA 
is not possible at this time because the reactions and their rates of 
pyrolysis flow have not been well studied. The only pyrolysis gas 
chemistry found in the open literature is the work by April et al. [4]. 
Because no other choice was available, in this work, we adopt the gas 
chemistry of April et al. [4], and we assume that the pyrolysis gas 
formed in the PICA char layer consists of nine gas species, 
specifically, H 2 , CH 4 , CO, C 2 H 6 , C 2 H 4 , C 2 H 2 , C0 2 , H 2 0, and C, 
which are the most abundant species for some chemical equilibrium 
conditions. The eight chemical reactions considered in this work and 
their rates are listed in Table 1. A detailed discussion of these 
reactions is presented in the paper by April et al. [4], 

k = AT~ s e ( - E l RT) (9) 

The specific heat of solid is input as a function of temperature for 
both virgin and fully charred material. In partially pyrolyzed zones 
(p c < p < p v ), the specific heat is obtained from the mixing rule 

c p = «' P v + f 1 “ T> c P c (!0) 

where 



The weighting variable r is the mass fraction of virgin material, in a 
hypothetical mixture of virgin material and char, which yields the 
correct local density. The thermal conductivity k is weighted in a 
similar manner as c p . 

III. Boundary Conditions 

Conditions at the ablating surface are determined by convective 
and radiative heating and by surface thermochemical interactions 
with boundary-layer gases. The surface energy balance equation is 
written in general convective transfer-coefficient form as follows: 

p e u e C H {H r - h w ) + m c (h c - h w ) + m g (h g - h w ) 

+ a w q rw - oeJJl - T%,) - r/ cw = 0 (12) 

The first term in Eq. (12) represents the sensible convective heat flux. 
The sum of the second and third terms in Eq. (12) is defined as the 
total chemical energy at the surface. The fourth and fifth terms are the 
radiative heat fluxes absorbed and reradiated by the wall. 



Time (s) 

Fig. 9 Three arcjet test conditions for group II. 
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respectively, and the last term, q cw , represents the rate of heat 
conduction into the TPS. If the surface is at the chemical equilibrium 
condition, ACE and MAT codes can be used to generate tables of 
normalized mass blowing rate, B' c (= '”/■ ). For a chemical 

nonequilibrium surface, the char recession rate m c can be computed 
through coupling FIAT and a computational fluid dynamics code with 
finite-rate surface chemistry capability, such as data-parallel line 
relaxation [17]. The enhanced version of FIAT code can read either 
the B' tables or the char recession rates, m c . 

A blowing correction accounts for the reduction in transfer 
coefficients due to the transpiration of gases front pyrolysis and 
surface ablation into the boundary layer. The blowing rate correction 
equation for convective heat transfer is 

C H ln(l + 2/LB') 


where A is the blowing reduction parameter, C H is the heat transfer 
coefficient for the ablating surface, and C H \ is the heat transfer 
coefficient for the nonablating surface. With A = 0.50, Eq. (13) 
reduces to the classical blowing correction for laminar flow [18]. 

IV. Test Cases 

The computations presented in this section focus on the in-depth 
thermal response of PICA heat-shield material. The first group of 
simulations is based on a lunar return trajectory for Orion, and the 
second group uses a representative variety of stream conditions for 
Crew Exploration Vehicle (CEV) TPS ADP arcjet testing conducted 
at NASA Ames Research Center. The effects of gas momentum 
transfer and pyrolysis gas homogeneous reactions on PICA in-depth 
thermal response are examined. In each group, we first solve the 
Darcy-Forchheimer momentum equation under the assumption of 





Fig. 10 Surface and in-depth temperatures predicted by baseline FIAT and by enhanced FIAT with chemical equilibrium chemistry for a) heat flux of 
900 W /cm 2 and pressure of 64 kPa, b) heat flux of 600 W /cm 2 and pressure of 26 kPa, and c) heat flux of 200 W /cm 2 and pressure of 7 kPa. 
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chemical equilibrium to understand how the inclusion of gas 
momentum transfer and gas kinetic energy affects the in-depth 
temperature prediction. Then, the gas momentum equation is solved 
with finite-rate homogeneous reactions, developed by April et al. [4] 
for pyrolysis gas flow, to study how the gas phase chemistry 
influences the PICA material thermal response. For each group, the 
predictions of PICA material thermal response using nine-species 
finite-rate gas chemistry, equilibrium gas chemistry, and frozen gas 
chemistry are presented and compared with those using the baseline 
FIAT/PICA model developed for CEV TPS ADP. 

A. Group I 

The aerothermal environment for the first group is shown in Fig. 1 . 
This flight environment is a dual convective heat pulse based on a 
proposed lunar return trajectory for Orion. The first peak convective 
heat flux is 245 W/cm 2 at 1 05 s, and the second peak is 85 W/cm 2 at 
523 s. The radiative heat flux is a single pulse with a maximum of 
204 W/cm 2 at 80 s. 


The gray lines shown in Fig. 2 are the PICA in-depth temperature 
histories at various depths (surface, 0.24, 0.51, 1.15, 1.76, 3.81, and 
5.08 cm) predicted by the regular version of FIAT using the baseline 
thermal response model, PICAv3.3, developed for CEV TPS ADP. 
The black lines are predictions of the enhanced version of FIAT, 
which solves the Darcy-Forchheimer momentum equation and the 
global mass conservation equation with, in this case, the assumption 
of chemical equilibrium for the pyrolysis gas flow. Both predictions 
are performed using a baseline 1 19-species pyrolysis gas model. In 
general, the enhanced version predicts slightly lower in-depth 
temperature than the regular version for locations inside the char and 
pyrolysis zones. At depths of 3.81 and 5.08 cm, at which the 
temperature is not sufficiently high for pyrolysis to be significant, the 
difference between these two predictions is very small. Consideration 
of the gas momentum equation does not appear to have significant 
impact on the thermal response prediction. The difference in 
temperature prediction is mainly driven by the introduction of kinetic 
energy of pyrolysis gas in the total energy balance equation [Eq. (1 )]. 




Fig. 11 Temperature and gas pressure profiles for a) heat flux of 900 W /cm 2 and pressure of 64 kPa at time equal to 30 s, b) heat flux of 600 W /cm 2 and 
pressure of 26 kPa at time equal to 50 s, and c) heat flux of 200 W/cm 2 and pressure of 7 kPa at time equal to 100 s. 
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The gas kinetic energy is not considered in the energy balance terms 
of the regular version of FIAT. One can conclude that the kinetic 
energy of pyrolysis gas is a relatively small component in the global 
energy balance of PICA material, and thus, the difference in predicted 
temperatures between two simulations is small. Thus, for flow 
environments studied in this case, ignoring the kinetic energy of 
pyrolysis gas should not have a significant impact on the accuracy of 
PICA material thermal response predictions. However, in a material 
thermal response model with finite-rate gas chemistry, an accurate 
prediction of the gas flow rate is required for the computation of 
chemical species concentrations. 

The pyrolysis gas mass injection rate vs time is presented in Fig. 3. 
The prediction from the regular version of FIAT (in gray) is just 
slightly higher than that from the enhanced version (in black). Instead 


of solving the gas momentum conservation equation and each 
individual gas species mass conservation equation, the regular 
version of FIAT code assumes that all pyrolysis gas flowing outward 
with neither resistance nor accumulation, and thus, the gas flow rate, 
can be obtained by simply solving the global gas mass conservation 
equation. Based on the results of the current analysis, this assumption 
is proved to be reasonable for the flight conditions considered here. 

To examine the effect of variation of permeability on PICA in- 
depth thermal response, computations with various values of perme- 
ability are performed. The results presented are for the permeabilities 
equal to baseline value (K), twice baseline (2 K), and half of baseline 
(0.5 A"). Figures 4a and 4b show the in-depth temperature and gas 
pressure profiles at 80 s, the first peak of surface temperature, and at 
515 s, the second peak of surface temperature, respectively. The 




e) 

Fig. 12 Gas density and solid density profiles for a) heat flux of 900 W /cm 2 and pressure of 64 kPa at time equal to 30 s, b) heat flux of 600 W/em 2 and 
pressure of 26 kPa at time equal to 50 s, and c) heat flux of 200 W/cm 2 and pressure of 7 kPa at time equal to 100 s. 
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temperature distributions are almost unaffected by the variation of 
permeability. The maximum temperature difference due to the 
variation of permeability is 20 K, which is too small to be seen in this 
plot. Thus, only the temperature distribution for baseline perme- 
ability is shown in this chart. The increase in pyrolysis gas pressure is 
about 15 to 20%, as the permeability is reduced by 50%, and 
conversely, the pyrolysis gas pressure decreases if the permeability 
increases. The predicted pyrolysis gas density and solid density are 
plotted in Figs. 5a and 5b for the same time steps (80 and 515 s). The 
effect of variation of permeability on gas density is similar to that on 
gas pressure. In Figs. 6a and 6b, the gas kinetic energy density and 
pyrolysis mass flow rate at time equal to 80 and 5 15 s are presented. 
As expected, flow speed and kinetic energy increase as the value of 


permeability increases. The total pyrolysis mass flow rate, which 
mostly depends on the resin decomposition rate, is not sensitive to the 
variation of permeability. 

The selection of chemical species is limited by the availability of an 
applicable finite-rate gas chemistry model. To better estimate the 
uncertainty associated with using a simple nine-species gas 
chemistry for predicting PICA material in-depth thermal response, a 
corresponding nine-species chemical equilibrium gas chemistry 
model was used. Both the equilibrium and the finite-rate models 
contain exactly the same nine chemical species. In Fig. 7a, the 
comparison of FIAT predictions between the baseline 119-species 
equilibrium chemistry model (in black), which is used in the regular 
version of the FIAT/PICA model, and the simple nine-species 
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Fig. 13 Kinetic energy density and mass flow rate profiles for a) heat flux of 900 W /cm 2 and pressure of 64 kPa at time equal to 30 s, b) heat flux of 
600 W /cm 2 and pressure of 26 kPa at time equal to 50 s, and c) heat flux of 200 W /cm 2 and pressure of 7 kPa at time equal to 100 s. 
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equilibrium chemistry model (in gray) is presented. The difference in 
predicted in-depth temperature profiles between the nine-species and 
119-species chemical equilibrium models is fairly small, except at 
depths of 1.15 and 1.76 cm, where the nine-species predictions are 
lower than those of 119-species predictions. This discrepancy is 
mainly due to differences in the gas enthalpy inside the pyrolysis 
zone. The difference in predicted in-depth temperature between these 
two equilibrium chemistry models is considered to be acceptable for 
the purpose of this parametric study. Since there is no other finite-rate 
gas chemistry model available, we have to accept the errors that may 
be generated as the result of using a simple nine-species model for the 
parametric studies performed in this section. The number of gas 
phase chemical species should be increased in a future study, when a 
more robust finite-rate homogeneous reaction model becomes 
available. 



Figure 7b presents the in-depth temperature profiles for the finite- 
rate (black lines), frozen (dark gray lines), and equilibrium (light gray 
lines) models with nine chemical species. The frozen chemistry 
solutions are obtained by turning off all the reactions in the finite-rate 
model. The greatest difference among these three models occurs at a 
depth of 1.76 cm. Generally speaking, the frozen chemistry has the 
highest temperature predictions, and the equilibrium chemistry has 
the lowest. This discrepancy is certainly due to the difference in the 
gas phase chemical reaction rates in each model. 

A comparison of temperature profiles between the nine-species 
finite-rate chemistry model (black lines) and the regular FIAT model 
(gray lines) is presented in Fig. 7c. The overall difference between 
these two predictions is not insignificant. At depths of 0.24 and 
0.51 cm, where PICA is fully charred most of the time, the regular 
FIAT model predicts a slightly lower temperature as compared with 
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Fig. 14 Surface and in-depth temperatures predicted by finite-rate, equilibrium, and frozen gas chemistry models for a) heat flux of 900 W/cm 2 and 
pressure of 64 kPa, b) heat flux of 600 W /cm 2 and pressure of 26 kPa, and c) heat flux of 200 W/cm 2 and pressure of 7 kPa. 
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the nine-species finite-rate FIAT, due to the difference in accounting 
for gas kinetic energy. At deeper locations, including 1.76, 3.81, and 
5.08 cm, the kinetic energy term is less important. The gas enthalpy is 
the main difference between these two models, and thus, the 
temperature predictions of the nine-species finite-rate model are 
slightly higher than those of the regular version. 

Figures 8a and 8b compare the mass concentrations for C0 2 , 
C 2 H 2 , CH 4 , H 2 , and CO predicted by the finite-rate (solid lines) and 
equilibrium (dashed lines) gas chemistry models at 80 and 515 s, 
respectively. The mass fractions for the other four species are not 
presented because of their relatively low concentrations. The 
pyrolysis gas flow is first produced inside the pyrolysis zone, expands 
into the char zone, reaches the front face, and then mixes with the flow 
surrounding the TPS system. For a chemical equilibrium model, the 
local chemical species concentrations depend on local temperature 



and pressure. However, for a finite-rate model, the individual 
chemical species concentration is driven by the production rate of 
each individual species, which is a time-dependent process. The 
history of local species concentrations, temperature, pressure, and 
flow speed all affect the extent to which the finite-rate chemistiy 
model approaches a chemical equilibrium condition. At time equal to 
80 s, the finite-rate model prediction gradually approaches chemical 
equilibrium as the pyrolysis gas flow moves near the front surface. 
Nonetheless, in this same region, at time equal to 515 s, the species 
concentrations predicted using the finite-rate model still have not 
reached those predicted using the chemical equilibrium model, 
because of relatively low gas density and pressure. We again note that 
this finite-rate model excludes heterogeneous reactions that may be 
important and that could increase the rate at which pyrolysis gas 
chemistiy approaches equilibrium conditions. 




Fig. 15 Mass fraction profiles predicted by finite-rate and chemical equilibrium chemistry models for a) heat flux of 900 W/em 2 and pressure of 64 kPa 
at time equal to 30 s, b) heat flux of 600 W/cm 2 and pressure of 26 kPa at time equal to 50 s, and c) heat flux of 200 W/cm 2 and pressure of 7 kPaattime 
equal to 100 s. 
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B. Group II 

Figure 9 shows three representative arcjet heating environments 
for the second group of FIAT simulations. The first case has a surface 
heat flux of 900 W / cm 2 , pressure of 64 kPa, and exposure duration 
of 30 s. The second case has a surface heat flux of 600 W/cm 2 , 
pressure of 26 kPa, and exposure duration of 50 s. The third case has a 
surface heat flux of 200 W/cm 2 . pressure of 7 kPa. and exposure 
duration of 100 s. The computations performed in this group are 
similar to those in group I. 

Figures 10a to 10c present temperature histories at various depths 
(surface, 0.24, 0.51, 1.15, 1.76, 3.81, and 5.08 cm) computed by the 
regular version of FIAT (gray lines) and the enhanced version of FIAT 
(black lines) using the baseline 119-species equilibrium chemistry 
model for all three arcjet simulation cases. The results again indicate 
that consideration of gas momentum transfer and gas kinetic energy 
only slightly decreases the predicted in-depth temperature. This result 
is obtained because the kinetic energy of pyrolysis gas is relatively 
small compared with other components of internal energy. 

The effect of permeability variations on gas pressure, density, and 
kinetic energy is also studied in this group. In Figs. 1 la to 11c, the gas 
temperature and pressure profiles at the end of each heat pulse for 
three permeability levels (0.5K, K, and 2 K) are presented. The gas 
density and solid density distributions at the end of the heat pulse are 
given in Figs. 12a to 12c, and the distributions of kinetic energy 
density and mass flow rate are shown in Figs. 13a to 13c. An increase 
in permeability reduces gas pressure and density; however, it 
increases gas flow speed and kinetic energy. These predictions 
confirm that gas mass flow rate and temperature are not particularly 
sensitive to the variation of permeability for PICA. 

A comparison of predicted temperature histories for the nine- 
species finite-rate (black lines), frozen (dark gray lines), and 
equilibrium (light gray lines) gas chemistry models is provided in 
Figs 14a to 14c for the three arcjet conditions. The observed 
discrepancy in the PICA in-depth temperature predictions among 
three chemistry models mainly results from the gas phase enthalpy. 
This is because the predicted chemical species concentrations 
strongly depend on the chemistry model implemented in the code, 
and the gas enthalpy is a function of chemical species concentrations. 
Inside the pyrolysis zone, the maximum difference in predicted in- 
depth temperature is 40 K. For this arcjet group, the effect of 
uncertainty on gas phase chemistry modeling of pyrolysis gas 
appears to have less influence on the overall PICA thermal response 
as compared with the first group of flight simulations, because of the 
relatively thin char layer. 

Figures 15a to 15c show the predicted mass fractions of COi, 
C 2 H 2 , CH 4 . H 2 , and CO for nine-species finite-rate (solid lines) and 
nine-species equilibrium (dashed lines) chemistry models at the end 
of each heat pulse. The trends are similar to what was seen in group 
I. For all three arcjet conditions, the results of the nine-species 
finite-rate chemistry model indicate that the pyrolysis gas flow is at 
or near chemical equilibrium in the region near the front surface 
of PICA. 


V. Conclusions 

A special version of fully implicit ablation and thermal response 
code was developed. The new capabilities implemented in this 
version of code include the time-dependent pyrolysis gas flow 
momentum equation with Darcy-Forchheimer terms and the pyrol- 
ysis gas species conservation equations with a nine-species finite-rate 
homogeneous chemical reaction model developed specifically for the 
phenolic impregnated carbon ablator. The total energy conservation 
equation was also enhanced to reflect these additions. Various 
parametric studies were performed for two groups of test cases. The 
first group used a heating environment for a proposed Orion lunar 
return vehicle, and the second group used three heating environments 
for arcjet tests conducted at NASA Ames Research Center. The 
central focus of these parametric studies was to understand the effect 
of pyrolysis gas chemistry on the in-depth thermal response. Thus 
finite-rate, equilibrium, and frozen homogeneous gas chemistry 
models for pyrolysis gas were examined. 


Results from the current parametric studies show that gas kinetic 
energy is a relatively small component in the internal energy balance 
equation, and neglecting gas kinetic energy has little negative impact 
on the accuracy of in-depth temperature predictions. The in-depth 
temperature predictions using a finite-rate or a frozen pyrolysis gas 
chemistry model are higher than those using a chemical equilibrium 
model. The results computed using a finite-rate gas chemistry model 
indicate that pyrolysis gas flow is at or near a chemical-equilibrium 
condition in the region near the front surface for all flow conditions 
studied in this work. Compared with the baseline model, this 
enhanced version of code tends to predict slightly lower temperatures 
in the char zone due to the effect of pyrolysis gas kinetic energy, and 
slightly higher temperatures in the pyrolysis zone as the result of 
finite-rate chemical reactions. However, overall, the differences in 
predicted in-depth temperatures are considered to be small. The 
presence of chemical nonequilibrium pyrolysis gas flow does not 
significantly alter the in-depth thermal response performance 
predicted using the chemical equilibrium gas model. Thus, the PICA/ 
FIAT material thermal response model developed by the Orion 
Thermal Protection System Advanced Development Project should 
be applicable for a fairly wide range of entry conditions with 
reasonably good accuracy. 
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